Is the diet cyclic phase-dependent in boreal vole populations?

Herbivorous rodents in boreal, alpine and arctic ecosystems are renowned for their multi-annual population cycles. Researchers have hypothesised that these cycles may result from herbivore-plant interactions in various ways. For instance, if the biomass of preferred food plants is reduced after a peak phase of a cycle, rodent diets can be expected to become dominated by less preferred food plants, leading the population to a crash. It could also be expected that the taxonomic diversity of rodent diets increases from the peak to the crash phase of a cycle. The present study is the first to use DNA metabarcoding to quantify the diets of two functionally important boreal rodent species (bank vole and tundra vole) to assess whether their diet changed systematically in the expected cyclic phase-dependent manner. We found the taxonomic diet spectrum broad in both vole species but with little interspecific overlap. There was no evidence of systematic shifts in diet diversity metrics between the phases of the population cycle in either species. While both species' diet composition changed moderately between cycle phases and seasons, these changes were small compared to other sources of diet variation-especially differences between individuals. Thus, the variation in diet that could be attributed to cyclic phases is marginal relative to the overall diet flexibility. Based on general consumer-resource theory, we suggest that the broad diets with little interspecific overlap render it unlikely that herbivore-plant interactions generate their synchronous population cycles. We propose that determining dietary niche width should be the first step in scientific inquiries about the role of herbivore-plant interactions in cyclic vole populations.

Gut microbiome biogeography in reindeer supersedes millennia of ecological and evolutionary separation.

Ruminants are dependent on their gut microbiomes for nutrient extraction from plant diets. However, knowledge about the composition, diversity, function, and spatial structure of gut microbiomes, especially in wild ruminants, is limited, largely because analysis has been restricted to faeces or the rumen. In two geographically separated reindeer subspecies, 16S rRNA gene amplicon sequencing revealed strong spatial structuring, and pronounced differences in microbial diversity of at least 33 phyla across the stomach, small intestine, and large intestine (including faeces). The main structural feature was the Bacteroidota to Firmicutes ratio, which declined from the stomach to the large intestine, likely reflecting functional adaptation. Metagenome shotgun sequencing also revealed highly significant structuring in the relative occurrence of carbohydrate-active enzymes (CAZymes). CAZymes were enriched in the rumen relative to the small and large intestines. Interestingly, taxonomic diversity was highest in the large intestine, suggesting an important and understudied role for this organ. Despite the two study populations being separated by an ocean and six millennia of evolutionary history, gut microbiome structuring was remarkably consistent. Our study suggests a strong selection for gut microbiome biogeography along the gastrointestinal tract in reindeer subspecies.

Faecal metabarcoding provides improved detection and taxonomic resolution for non-invasive monitoring of gastrointestinal nematode parasites in wild moose populations.

BACKGROUND: Although wild ungulate populations are heavily monitored throughout Europe, we understand little of how parasites affect population dynamics, and there is no systematic, long-term monitoring of parasite diversity and parasite loads. Such monitoring is in part hampered by a lack of time- and cost-effective assay methodologies with high sensitivity and good taxonomic resolution. DNA metabarcoding has been successfully used to characterize the parasitic nemabiome with high taxonomic resolution in a variety of wild and domestic hosts. However, in order to implement this technique in large-scale, potentially non-invasive monitoring of gastrointestinal parasitic nematodes (GIN), protocol optimization is required to maximize biodiversity detection, whilst maintaining time- and cost-effectiveness. METHODS: Faecal samples were collected from a wild moose population and GIN communities were characterized and quantified using both parasitological techniques (egg and larva counting) and DNA metabarcoding of the ITS2 region of rDNA. Three different isolation methods were compared that differed in the volume of starting material and cell lysis method. RESULTS: Similar nematode faunas were recovered from all samples using both parasitological and metabarcoding methods, and the approaches were largely congruent. However, metabarcoding assays showed better taxonomic resolution and slightly higher sensitivity than egg and larvae counts. The metabarcoding was not strictly quantitative, but the proportion of target nematode sequences recovered was correlated with the parasitologically determined parasite load. Species detection rates in the metabarcoding assays were maximized using a DNA isolation method that included mechanical cell disruption and maximized the starting material volume. CONCLUSIONS: DNA metabarcoding is a promising technique for the non-invasive, large-scale monitoring of parasitic GINs in wild ungulate populations, owing to its high taxonomic resolution, increased assay sensitivity, and time- and cost-effectiveness. Although metabarcoding is not a strictly quantitative method, it may nonetheless be possible to create a management- and conservation-relevant index for the host parasite load from this data. To optimize the detection rates and time- and cost-effectiveness of metabarcoding assays, we recommend choosing a DNA isolation method that involves mechanical cell disruption and maximizes the starting material volume.

Issues of under-representation in quantitative DNA metabarcoding weaken the inference about diet of the tundra vole Microtus oeconomus.

During the last decade, methods based on high-throughput sequencing such as DNA metabarcoding have opened up for a range of new questions in animal dietary studies. One of the major advantages of dietary metabarcoding resides in the potential to infer a quantitative relationship between sequence read proportions and biomass of ingested food. However, this relationship's robustness is highly dependent on the system under study, calling for case-specific assessments. Herbivorous small rodents often play important roles in the ecosystem, and the use of DNA metabarcoding for analyses of rodent diets is increasing. However, there has been no direct validation of the quantitative reliability of DNA metabarcoding for small rodents. Therefore, we used an experimental approach to assess the relationship between input plant biomass and sequence reads proportions from DNA metabarcoding in the tundra vole Microtus oeconomus. We found a weakly positive relationship between the number of high-throughput DNA sequences and the expected biomass proportions of food plants. The weak relationship was possibly caused by a systematic under-amplification of one of the three plant taxa fed. Generally, our results add to the growing evidence that case-specific validation studies are required to reliably make use of sequence read abundance as a proxy of relative food proportions in the diet.

Comparing three types of dietary samples for prey DNA decay in an insect generalist predator.

The rapidly growing field of molecular diet analysis is becoming increasingly popular among ecologists, especially when investigating methodologically challenging groups, such as invertebrate generalist predators. Prey DNA detection success is known to be affected by multiple factors; however, the type of dietary sample has rarely been considered. Here, we address this knowledge gap by comparing prey DNA detection success from three types of dietary samples. In a controlled feeding experiment, using the carabid beetle Pterostichus melanarius as a model predator, we collected regurgitates, faeces and whole consumers (including their gut contents) at different time points postfeeding. All dietary samples were analysed using multiplex PCR, targeting three different length DNA fragments (128, 332 and 612 bp). Our results show that both the type of dietary sample and the size of the DNA fragment contribute to a significant part of the variation found in the detectability of prey DNA. Specifically, we observed that in both regurgitates and whole consumers, prey DNA was detectable significantly longer for all fragment sizes than for faeces. Based on these observations, we conclude that prey DNA detected from regurgitates and whole consumers DNA extracts are comparable, whereas prey DNA detected from faeces, though still sufficiently reliable for ecological studies, will not be directly comparable to the former. Therefore, regurgitates and faeces constitute a useful, nonlethal source for dietary information that could be applied to field studies in situations when invertebrate predators should not be killed.